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ABSTRACT 

The time anomaly of neutrino events observed in the KARMEN detector 
at the ISIS facility is discussed. We show that beam-correlated neutron 
events are not completely suppressed by the KARMEN lower energy cut and 
could cause the observed deviation of the measured time distribution from 
an exponential curve with r = 2.2 fis. 
PACS: 13.20. Cz, 14.60.Pq, 14.80.-j 
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The KARMEN collaboration investigates reactions induced by neutrinos pro- 
duced in the pulsed neutron spallation source of the ISIS facility at the 
Rutherford Appleton Laboratory. In 1995, they reported ||1| an anomaly 
in the time distribution of their events clustered around 3.6 /is after beam- 
on-target with visible energy between 11 and 35MeV. Further evidence of 
events exceeding the expected exponential distribution characterized by the 
muon lifetime has since been reported ||2|. The speculative explanation was 
that these events could be generated by hitherto unobserved neutral massive 
X-particles originating from the ISIS target via the rare pion decay process, 

7r+ ^ 12+X. (1) 

From the kinematics of the two-body decay of a pion at rest, and the known 
masses of the pion and the muon 0, the mass of X was calculated to be 
33.9 MeV using the time-of-flight information and the mean flight-path of 
17.5 m from the ISIS target to the KARMEN detector. Since the detected 
energy of the anomalous events is found to be much greater than the kinetic 
energy of the X-particle, these events are postulated to originate from X- 
particles which decay with visible energy within the detector via the reaction 
X —>■ e+e^z/. The KARMEN collaboration estimated the branching fraction 
for the decay mode (1) as a function of the X-lifetime, with values of the 
branching fraction ranging over many decades down to ~10~^^ |[l|. 

Soon after the initial KARMEN publication various hypotheses were 
put forward as to the nature of the X-particle and several searches 
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undertaken to look for the rare decay (1) 10-0], all with negative results. 



More recently, additional experiments |[13|-[T5|| have set even more stringent 
limits for the branching fraction of the decay mode (1). 

In this letter we present arguments for the existence of a beam-correlated 
neutron background in the KARMEN detector that could contribute to an 
excess of event^. Our basis is that the KARMEN detector can respond to 
the 10.8 MeV 7-energy from neutron capture by ^"^N. Thermal neutrons will 
originate from moderation of fast neutrons (T„ < 10 MeV) in the liquid scin- 
tillator. These fast neutrons are produced from medium energy interactions 
(T„ > 100 MeV) in the KARMEN iron shield. 

The liquid scintillator (PPP) used in the KARMEN detector consists of 
75% (volume) paraffin oil, 25% pseudocumene and 2g/l PMP (1-phenyl- 



3-mesityl-2-pyrazoline) [|T^, |T7|. The chemical formulae of the components 
are (i) paraffin oil, taken as C11H24, (ii) pseudocumene C9H12, (iii) PMP 
C18H20N2. The nuclear densities in the liquid scintillator material (in units 
of nuclei/A^) are 0.071 for hydrogen, 0.036 for carbon and 9.1 • 10~^ for 
nitrogen leading to the ratio of atomic abundances of 

uh -.nc-.nN ^ 7800 : 4000 : 1. (2) 

It is this nitrogen that can capture neutrons via the reaction 

^^N(n,7)^^N (3) 

with an integral 7-energy of 10.8 MeV. The thermal neutron cross-section for 
^The number of excess events reported in Ref. pi is 83 ± 28. 



this reaction is a ~75 mb. The lower energy cut of events in the KARMEN 



analysis is llMeV with an energy resolution of cte/E = 11.5%/v-E 
{E in MeV), which is 3.5% or 380 keV for E = llMeV. Thus, the 10.8 MeV 
7-events from reaction (3) will be accepted with an efficiency of about 30%; 



we note that in the energy distribution of the excess events [0, about 40% 
of the anomalous events, which are selected from an energy window between 
11 and 35 MeV, are clustered at about 11 MeV. If such a beam-correlated 
neutron background is present, the subtraction of purely cosmic background 
events is not sufficient. 

The data fit published by the KARMEN collaboration gives a value for 
the time constant = (2.62 ± 0.18) /iS ||l| which differs by more than two 
standard deviations from the expected value of the muon lifetime (2.19703 ± 
0.00004 /is). We have refitted the KARMEN data of Ref. g with an addi- 
tional background of the form: 



y 



it) = A- e-'l^^ + 5, (4) 



and fixing to the value of 2.2 /xs. The result of our fit is shown in Fig. 1. 
This fit has a of 16.4 for 17 DOF and the values of the two free parameters 
found are A = (271 ± 16) events/0.5/is, B = (7.2 ± 3.0) events/0.5/is. The 
value for B, greater than zero by 2.4 standard deviations, can be interpreted 
as an indication for the presence of an additional background with 137 ± 57 
event s0. 

2 A similar fit was made to the combined data from KARMENl and KARMEN2 ITll. 
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From the thermal capture cross-sections of hydrogen (333 mb), carbon 
(3.5 mb) and nitrogen (75 mb for n,7 and 1810 mb for n,p) and using the 
ratio of atomic abundances (2), we deduce the capture cross-section of the 
scintillator to be about 220 mb and the probability for reaction (3) to be 
P ^ 3 ■ 10"^. This means that 1/P ~ 3 ■ 10^ thermal neutrons would be 
needed to create photons of reaction (3) with 10.8 MeV of deposited visible 
energy. In order to pass the 11 MeV cut in the KARMEN analysis, this 
number has to be a factor of 3 larger, i.e. about 1 ■ 10^. 

Because of the neutron absorbers around each scintillator element (gadolin- 
ium and borated polyethylene) ||l|, thermal neutrons cannot enter the KAR- 
MEN detector and hence have to be produced by moderating fast neutrons 
in the scintillator fluid itself. Since the hydrogen density in the PPP is sim- 
ilar to that for water, the moderation characteristics should be alike. 'Age 
theory' ||19[ predicts an average time for thermalization of a few microsec- 
onds. The time variation of the background signal will be complicated by 
the 1/v-tailback of the capture cross-section. Thus, the contributions to the 
background signal, will start well before thermalization is complete. 

The events of Fig. 1 originate from 6560 C, i.e. 4 ■ 10^^ protons hitting the 
ISIS target ||l|. Each proton produces about 0.05 medium energy neutrons 



per steradian at 90 degree production angle ||20|. The KARMEN detector 
covers 5 • 10^^ of the solid angle. The material between the target and the 
KARMEN detector provides about 4600 g-cm'^ of shielding ||I],|16|,|2g,|21|, and 
the shielding length for medium energy (above ~ 100 MeV) neutrons is in the 



range 120 to 160g-cm~^ [0]. This leads to the estimate that between 5 ■ 10'' 
and 3 ■ 10'^ medium energy neutrons enter the KARMEN detector. Medium 
energy neutrons contribute 0.4% of the total neutron flux escaping from an 
iron shield ||2^, and 98 % of all neutrons have an energy below 10 MeV. With 
these numbers, we get the contributions to our probability estimate of a 
beam-correlated neutron background to be detected in the KARMEN time 
window (10/is) listed in Table 1. 

From this table, we calculate the suppression factor of beam-correlated 
neutrons to be between I-IO"^*^ and 2-10~^° for the two extremes of the shield- 
ing length. Thus, for a total number of protons (4 ■ 10^^), the corresponding 
expected number of neutron events originating from a beam-correlated back- 
ground is between 0.06 and 800 which brackets the number of background 
events (137 ± 57). 

For a more reflned background estimate, one has to consider that in a 
pulsed neutron spallation source, a beam-correlated neutron background has 
a time structure originating from the duty cycle of the primary proton beam 
and the kinetic energy of the neutrons on their way through the shielding. 
The neutron capture rate of gadolinium absorber dissolved in water was 
measured as a function of time after injection of a fast neutron pulse ||23|| and 
corresponds to a linear rise from zero to about 6 fis followed by a constant rate 
at longer timesR. The time variation of the neutron capture rate in the liquid 



^In fact, an exponential fall with a time constant of ~ 100 |16| (the thermal neutron 
lifetime in the KARMEN scintillator) is to be expected. 
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scintillator is expected to be similar. This behavior has been approximated 
by the function 

y{t) = A ■ e(-*/2-2M«) + B ■ t/r, (5) 

for t<Tc and 

yit) = A ■ e(-*/2-2M^) + B (6) 

for t>Tc. Here, A is the initial rate of the exponential decay, B is the beam- 
correlated background from reaction (3), and Tc is the linear risetime of the 
neutron capture rate. A fit to the data of Ref. results in a of 15.5 
for 16 DOF. The fitted values of the free parameters were A = (277 ± 14) 
events/0. 5/is, B = (7.3 ± 2.8) events/0. 5/is, and Tc = (3.4 ± 0.5) /xs. This 
fit, shown in Fig. 2, leads to 120 ± 42 neutron events from a beam-correlated 
background. Although the statistical significance is similar to the previous 
fit, it gives a possible explanation for an enhancement in the spectrum due 
to the neutron capture rate in the scintillator reaching its maximum some 3 
to 4/is after beam on target. 

After the upgrade of the KARMEN shielding (1996/97), any beam-correlated 
background from neutrons should be further suppressed due to additional 
polyethylene, veto scintillation counters, and iron shielding. Thus, the cor- 
responding 'time anomaly' in the KARMEN data (|,|18| before 1997, should 
be less significant in the data after 1997, as was presented in Ref. p^ . 

In conclusion, it can at present not be ruled out that the number of 
KARMEN excess events and their time distribution according to Ref. ||l| are 
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consistent with the hypothesis of originating from a beam-correlated neutron 
background involving neutron capture in nitrogen via reaction (3). A more 
elaborate analysis of the time structure of such a neutron background 



is beyond the scope of this letter, however, we would like to mention the 
existence of deep minima in the total cross section of iron, e.g. around 81 keV 
p5| . At this energy, neutrons travel through iron with a mean free path of 
several meters with a velocity of 4m//is and have a TOF of ~ 2 fis from 
the surface of the ISIS biological shielding to the KARMEN detector. These 
neutrons have a time structure similar to that of the proton beam on the ISIS 
target and so could lead to a localized enhancement of the beam-correlated 
background event-rate. 

The authors thank G. H. Eaton, R. Eichler, R. Frosch, and V. Markushin 
for valuable discussions and interest. This work was supported by the Schweiz- 
erischer Nationalfonds. 
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Figure Captions 

Figure 1: 

KARMEN data of Ref. g fitted with the function y(t) = A ■ e-*/2-2M« + B. 
The smooth fine represents the fitting function. The dashed fine is an expo- 
nential with = 2.2 /xs, and the dotted line represents an additional back- 
ground. The values of the fitted parameters are A = (271 ± 16) events/0. 5/is, 
B = (7.2 ± 3.0) events/0.5/is; the of the fit is 16.4 for 17 DOF. 

Figure 2: 

KARMEN data of Ref. g fitted with the function y(t) = A-e'^/^-^f"' + B -t/r^ 
for t < Tc and y(t) = A - e~*/2-2Ms _|_ q fQ^. t > Tc- The smooth line represents 
the fitting function. The dashed line is an exponential with = 2.2 /is, and 
the dotted line represents the time dependence of our estimate of neutron 
capture in the KARMEN scintillator from reaction (3). The values of the 
fitted parameters are A = (277 ± 14) events//is, B = (7.3 ± 2.8) events//is, 
and Te = (3.4 ± 0.5) /is; the of the fit is 15.5 for 16 DOF. 
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Figure 1: 




time rusl 
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Figure 2: 




time Tasl 
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Table 1: Factors used for estimating the beam-correlated neutron background 



Number of medium energy neutrons directed 
towards KARMEN 


0.05/sr/proton 


Solid angle for the KARMEN detector 


0.0628 sr 


Shielding factors for 4600 g/cm^ and shielding lengths 
A of 120g/cm^ — 160g/cm^ 


2- 10-1^-3 - 10-^3 


Fraction of neutrons escaping the KARMEN steel shield 
in the medium energy region (0.4%), 
in the fast neutron region (98 %) ; the number of fast 
neutrons per medium energy neutron is 


0.98/0.004 ^ 250 


Shielding effect for fast neutrons of the borated 
polyethylene layer and the active outer shield 


0.25 


Fraction of neutrons captured by PPP 


0.78 


Probability of the capture channel ^^N(n, 7)^^N 


3 • 10-^ 


Detection efficiency for 10.8 MeV (energy cut llMeV) 


0.3 


Fraction of captures within the KARMEN time window 


0.05 
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